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Abstract. It is proven that, under mild physical assumptions, an isolated 
stationary relativistic perfect fluid consists of a finite number of cells fibred 
by invariant annuli or invariant tori. For axially symmetric circular flows it is 
shown that the fluid consists of cells fibred by rigidly rotating annuli or tori. 



1. Introduction 

The purpose of this paper is to extend some results in the non-relativistic the- 
ory of Euler flows to the theory of stationary perfect fluids in general relativity. 
The latter play an important role in understanding the Einstein field equations in 
the presence of matter, while the former sits within the larger field of topological 
hydrodynamics. In particular, our main result is: 

Theorem 1: Consider an analytic stationary perfect fluid flow on an analytic 
spacetime ( A M , A g ) with proper pressure p and density fi satisfying Einstein's field 
equations, which has the following attributes: 

(a) it is spatially compact, 

(b) it satisfies the weak energy condition: p + fj, > 0, 

(c) it is isentropic: [i = /i(p), 

(d) it is non- Beltrami (see Definition 3 below). 

Then the spatial fluid region is divided into a finite number of cells of two types: 
if the cell does not intersect the boundary it is fibred by tori invariant under the 
flow, if the cell intersects the boundary it is fibred by annuli invariant under the 
flow. The flow lines on a torus are either all closed or all dense, while the flow 
lines on an annulus are all closed. 
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To prove this theorem we find a pair of commuting vector fields that generate the 
invariant surfaces. These invariant surfaces have less geometric significance than 
the surfaces of transitivity of axial symmetry, but the theorem states that they are 
(generically) topologically restricted. No more than the conservation equations of 
general relativity are needed to prove this. 

We apply this theorem in the case where the fluid is axially symmetric and the 
fluid 4-velocity lies in the plane spanned by two Killing vectors, so called circular 
flows: 

Theorem 2: A circular fluid flow is made up of cells foliated by rigidly rotating 
annuli or tori. 

In the next section we summarise the background material of stationary perfect 
fluids needed to prove the Structure Theorem. This includes a reformulation of the 
conservation equations for a relativistic stationary perfect fluid along the lines of 
the Euler equations for a steady incompressible fluid flow. Section 3 provides the 
definition of relativistic Beltrami flows and the proof of Theorem 1 while section 4 
applies the theorem to circular fluid flows. The final section contains a discussion 
of the results and their implications. 

2. Stationary Perfect Fluid Flows 

A lorentz 4-manifold ( 4 M, 4 #) is said to be stationary if it admits a timelike 
Killing vector We can choose coordinates (x° — t, x l ) on a stationary spacetime 
so that = 5q and the metric takes the i-independent form 

ds 2 = g io dx l dx j - e u {dt + Adx 1 ) 2 . 

Here and throughout Latin indices run through 1,2,3, while Greek indices run 
through 0,1,2,3. 

Locally we can quotient out by the action of the Killing vector on 4 M to get a 
riemannian 3- manifold (M, g) , which we think of as the spatial slice of the space- 
time. We use this riemannian metric to raise and lower all spatial tensor indices. 
However, for tensors defined on the whole spacetime we will continue to use 4 g. For 
example, since is a 4- vector 

Ui = 4 gi „U» = 4 9ij W + 4 g t0 U° = 9iJ U* - e"A t A 3 W + 4 g l0 U°, 

while for a 3- vector V 1 , Vj = g%jV^ . 

For later use, we note that for any f-independent 4- vector the divergence with 
respect to the lorentz and riemannian metrics are related by the simple formula 

4 V^ = V l W + \u),ifP. 
For any 4-vector we note the following useful relations 

U° = -e^Uo - A l U t U = -e^(U Q + A l U t ) 

U l = g ij Uj - A l U Q U, = 9ij U j + AiU , 

1 we have that 

er^U 2 = l + g ij U i U j . 



and if U U U^ = - 
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Consider a stationary spacetime ( 4 M, 4 g) on which we have a perfect fluid with 
4- velocity U v , pressure p, density /i and energy momentum tensor 

T a p = (p + n)U a Ufj + p i g aj3 , 

and U V U V = -1. 

The assumption that the flow is steady on a stationary spacetime means that 
U v , p and /i are all independent of t. 

For such a fluid, the conservation equations 4 V A T^ M = read 

4 V„p + 2(p + v)U a 4 V [a U v] + (U a 4 V Q (p + M ) + (p + M) 4 V q [/ q ) £7„ = 0. 

(2-1) 

Multiplying this across by U v we see that 

(p + ^) 4 V a U a + U a4 V a ^^0. (2.2) 

Substituting this back into (2.1) we get that 

4 W uP + 2(p + m)C/ q 4 V [Q [/, y] + 4 V Q p = 0. 

Now taking the v = component we get 

C/of'ViP + (p + tilTViUa = 0, (2.3) 

while taking the v = i component we get 

V tP + 2(p + ^tPVyC^ - (p + tiU*ViU Q + UiWVjP = 0. 

Multiplying this across by ?7o, using (2.3) and the fact that Z7o^° = — 1 — t/iC/ 2 we 
find that 

U V lP +( P + m)V«C/ + 2(p + riUStfVuiUo 1 ^) = °- ( 2 - 4 ) 

Equations (2.2) and (2.4) are the relativistic Euler equations governing a steady 
fluid perfect fluid flow on a stationary spacetime. 

Definition 1. A perfect fluid is isentropic if it has an equation of state /i = nip). 

In what follows we will assume that the perfect fluid is isentropic. If we introduce 
the real-valued functions a and j3 on M given by 



> + fJ, 

then the conservation equations can be written 

V l U l + ^U t V l uj + WV l \n(f3) =0 (2.5) 

V, ln[aU ] + 2U WV y {U^U t] ) = 0. (2.6) 

Definition 2. An isolated stationary body is one for which the spatial set V C M 
on which the energy momentum tensor is non-zero is compact, and the Darmois 
conditions hold on the boundary [1]. 

For a perfect fluid (with p + fi ^ 0) being isolated implies the vanishing of the 
pressure and tangency of the fluid flow at the boundary. 
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3. The Structure Theorem 
Consider the two vector fields 

X 1 = U a W Y l = Uo/3- 1 e-^r 1 ijk V u (U^U k] ) , 

where w is the hodge star operator: 

^ijk _ \ c ijk 

e ijk bgjng the 3-dimensional Levi-Civita permutation symbol. 

Definition 3. We say that the flow is Beltrami if Y % = XX 1 for some real- valued 
function Aon K 

The reason for such a definition becomes clear from the following structure the- 
orem for rclativistic iscntropic flows: 

Theorem 1. Consider a stationary analytic perfect fluid flow on an analytic space- 
time that: 

(a) is isolated, 

(b) satisfies the weak energy condition p + /i > 0, 

(c) is isentropic: ^ = fi(p), 

(d) is non-Beltrami. 

Then the spatial fluid region is divided into a finite number of cells of two types: 
if the cell does not intersect the boundary it is fibred by tori invariant under the 
flow, if the cell intersects the boundary it is fibred by annuli invariant under the 
flow. The flow lines on a torus are either all closed or all dense, while the flow 
lines on an annulus are all closed. 

Proof. The key to the proof is to show that the vector fields X 1 and Y 1 commute. 
To show this we use the 3-dimcnsional vector notation 

U x V = n^UyV^ curl(ff) = V^ b U k] ^ 



dW(U) = V i U i gra d(p)= ff «gA. 
The Euler equations (2.5) and (2.6) say that 

div(X) = X ■ grad[-iw + ln(C/ ) - m(/?)] 

grad {ln[aU }) -Xx curl(f) = 0, 

where we have introduced the 3- vectors W l = U l and T l = g lJ U Q ~ 1 Uj = U ~ 1 W l + 
A 1 . To compute the commutator of X 1 and Y % we utilise the vector identity 

[X,Y] = curl(f x X) + Xdiv(Y) - Ydiv(X). 

With the help of the Euler equations and the vector identities curl(grad /) = and 
div(curl(T)) = 0, this says that 

[X,curl(f)] = curl(curl(f) x X) + Xdiv(curl(f )) - curl(f )div(X) 
= -curl(f )X ■ grad[-iw + ln(U ) - ln([3)]. 
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This can be rearranged to 

[X,Y] = [X,U o 0- 1 e-i u cuil(f)] = 0, 

as required. 

The rest of the proof now follows that of the non-relativistic case and therefore 
we only sketch it below - details can be found in [2] . 

By analyticity, the critical set K of the function \n(aU ) forms an semianalytic 
subset and M\K has a finite number of cells. 

Assume that X % and Y l are not everywhere collinear. We can conclude that at 
non-critical values for ln(af7o), the level sets form an oriented smooth 2-manifold 
with an R 2 action. If the level set does not intersect dV, it forms a torus with all 
flow lines closed or all flow lines dense. If it does intersect the boundary, it must 
form an annulus with closed flowlines. 

□ 

Note: A similar Structure Theorem holds in the case of dust, where the pressure 
is identically zero. 

4. Axially Symmetric Flows 

When one considers the definition of axial symmetry for a stationary field, there 
are a number of reasonable possibilities [3] [4] [5]. For our purposes we use the 
following definition: 

Definition 4. A stationary space-time ( 4 M, 4 <?), with timelike Killing vector is 
axially symmetric if there exists a spacelike Killing vector rf that has closed orbits 
and commutes with 

Now suppose we have a stationary, axisymmetric spacetime ( 4 M, 4 g) as defined 
above. The simplest type of flow appears to be one where the spatial fluid flow U l 
is parallel to rf . 

Definition 5. A fluid flow on a stationary axisymmetric spacetime is circular if 
the fluid 4- velocity U v lies in the plane spanned by the Killing vectors and rf . 
In addition we assume that the spacetime has a regular axis of symmetry [6] [7] . 

Proposition 1. A circular flow is Beltrami iff X = 0. 

Proof. Working locally there exist coordinates (x a , <f), t) such that the Killing vectors 
are 

-._d_ ^_d_ 

We now utilise a remarkable theorem [6] [7] that says that for a circular fluid flow, 
the 2-planes orthogonal to the Killing orbits generated by (f and rf are integrable. 
Thus the metric can be reduce to 

ds 2 = h ab dx a dx b + H 2 d(j) 2 - e"{dt + A d<f)) 2 , (4.1) 

where all of the components of the metric are independent of both t and (f>. 
Circularity means that the fluid flow vector is given by 

U" = U° + £V) , (4.2) 

for some function f2. In particular, the only non- vanishing component of the spatial 
vector X is 

X 4, = U U°fl. 
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On the other hand 

Y* = tfo/rV^detfo))"^ 6 ^ (U^W C] + A c] ) = 0. 
This means that X and Y are orthogonal, and so the flow is Beltrami iff Y — 0. □ 

Virtually all results on stationary axially symmetric perfect fluids, whether exact 
solutions, asymptotic behaviour or a priori estimates, pertain only to circular flows. 
Thus the metric has the form (4.1) and the fluid velocity (4.2). 

Definition 6. The function defined in equation (4.2) is called the angular velocity 
of the fluid. A spinning fluid body with constant angular velocity is said to be in 
rigid rotation. If the angular velocity is not constant, the fluid is said to have 
differential rotation. 

Theorem 2. A circular flow consists of cells foliated by rigidly rotating annuli or 
tori. 

Proof. First, we gather together some of the relationships enjoyed by the angular 
velocity, the spatial fluid velocity vector W l , the 4-components of the fluid U° and 
Uq, and the metric functions H, e u and A. 



U a = -(l + Atye^U = y/l + \W 

\W\ 2 



= u°n Wa = 



{u°n) 2 u i + An 

A=- o+TTnTTTrl (4-3) 



n u a u n 

where we have introduced \W\ 2 = gijW l W J ' . 

Now, introduce coordinates (p,z) on the 2-space orthogonal to the group orbits, 
so that the commuting vector fields are 

Here p labels the invariant surfaces and (z,(j>) are coordinates on the invariant 
surfaces, analogous to the action angles of classical mechanics. In these coordinates 
the Euler equations dramatically simplify to 

The first two of these, along with the identity (4.3), tell us that the angular velocity 
fi is constant on each invariant surface. Thus each annulus or torus rotates rigidly, 
although there may be differential rotation between the surfaces. In addition, U°Uo 
is constant on each invariant surface. □ 
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5. Discussion 

In topological (non-relativistic) hydrodynamics there is an equivalent Structure 
Theorem for steady incompressible perfect fluid flows on a riemannian 3-manifold 
[2]. The dynamics of the fluid is described by four equations: the three Euler 
equations and the incompressibility condition. In this setting a fluid is Beltrami 
if the velocity vector is everywhere collinear with it's curl. In contrast to the 
regular behaviour of non-Beltrami fluids guaranteed by the Structure Theorem, the 
flowlines of Beltrami fluids can be arbitrarily knotted [8] [9] [10] [11]. Generically a 
fluid will be non-Beltrami, but Beltrami fluids do arise, for example, as "force-free" 
fields in magneto-hydrodynamics [12]. 

In the relativistic setting, stationary perfect fluids have long been studied and 
while much progress has been made, a convincing overview has yet to emerge. 
The literature in this area is rather extensive and we mention just some topics 
pertinent to our result - for a review see [13]. In particular, there has been much 
progress in finding and investigating local solutions that satisfy various simplifying 
assumptions. These can range from the weak (like energy conditions) to the strong 
(like cylindrical symmetry), or a combination of both. 

Without doubt the most commonly imposed conditions are those of axial sym- 
metry and circular flow, usually supplemented by some other restriction [14] [15] 
[16] [17] [18]. This setting comes closest to Newtonian theory and thus we are 
entitled to attempt to physically interpret the geometric quantities that arise. As 
has been shown, however, even in this case, care must be taken [19]. A further 
down side is that it is not clear which of the properties these local solutions possess 
are just products of their specialty. Thus, failure to match globally to the vacuum 
Kerr solution, a difficulty that every known stationary perfect fluid possesses, may 
simply be an artifact of their axial symmetry or circularity. 

Aside from exact solutions, there has been progress on general matching condi- 
tions with the Kerr solution [20] [21], as well as results on a priori estimates on the 
equatorial radius [22] , positivity of angular momentum density [23] and asymptotics 
[24]. In all cases axial symmetry is assumed, or the fluid is assumed to have zero 
pressure, i.e. the body consists of dust. 

It is within this context that our Structure Theorem above should be viewed. 
For here, the assumptions made on the fluid arc nothing more than one would 
impose on any physically reasonable model for an isolated stationary fluid body. 
Moreover, such a stationary perfect fluid is precisely what is required to model the 
equilibrium states of rotating fluid bodies [25] [26] . 

Examples of both types of cells in the Theorem have appeared in the literature. 
The remarkable perfect fluid source for the NUT metric discovered by Lukacs [27] 
has flowlines given by the Hopf fibering of S 3 , which lie on invariant tori. The 
second type of cells arise, for example, when one considers cylindrically symmetric 
flows [28] [29] [30] [31]. The decomposition of the full Einstein field equations 
parallel and transverse to the foliation is worthy of careful consideration. 

It is hoped that the Structure Theorem can yield a new approach to finding a 
source for the Kerr metric. It indicates that the simplest such source must be a solid 
torus and the fact that the singularity in the maximally extended Kerr solution is 
ring shaped supports such a conclusion. 

Another point worthy of further consideration is the Beltrami condition. In 
the non-relativistic case, Beltrami flows exhibit chaotic dynamics and are closely 
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related to contact geometry. The natural definition in the relativistic case afforded 
us by the Structure Theorem is a good deal more complex and its relationship 
with contact structures is unclear. It nevertheless raises the question as to whether 
relativistic effects "tame" Beltrami flows, or whether chaotic regimes persist. 
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